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Using 1310.6 x 10° J/2) and 447.9 x 10° (3686) events collected with the BESIII detector at the 
BEPCII efe collider, the branching fractions and the angular distributions of J/w and w(3686) 
decays to AA and X?” final states are measured. The branching fractions are determined, with 
much improved precision, to be 19.43 + 0.03 + 0.33, 11.64 + 0.04 + 0.23, 3.97 + 0.02 + 0.12 and 
2.44 + 0.03 + 0.11 for J/y — AA, J/y > XPE’, y(3686) — AA and (3686) — X°S°, respectively. 
The polar angular distributions of 7(3686) decays are measured for the first time, while those of J/w 
decays are measured with much improved precision. In addition, the ratios of branching fractions 


B(a)(3686) AA) and B(wb (3686) +D° 5°) 
B(J/b—AQ) B(J/p— 5050) 


are determined to test the “12% rule”. 


PACS numbers: 12.38.Qk, 13.25.Gv, 23.20.En 


I. INTRODUCTION throughout the text), take place through annihila- 

tion of the constituent cc quark pair into either a vir- 

Two-body baryonic decays of ọ mesons (w de- tual photon or three gluons, and they provide a good 
notes both the J/i and p (3686) charmonium states laboratory for testing Quantum Chromodynamics 


201711.02285v1 


chinaXiv 


(QCD) in the perturbative energy regime and study- 
ing the properties of baryons {a}. Perturbative QCD 
(pQCD) predicts that the ratio of branching frac- 
tions between the J/w and w(3686) decaying into a 
given hadronic final states follows the “12% rule” B 


= Besoin- a (12.4 + 0.4)%. 
Byjpsiti- 
(1) 


The violation of this rule was first observed in the 
decay of w into the final state px, which is well 
known as the “pz puzzle” Ep and the rule has been 
subsequently further tested in a wide variety of ex- 
perimental measurements. Reviews of the theoret- 
ical and experimental results [i] conclude that the 
current theoretical understanding, especially for the 
w decays into baryon-antibaryon pair final states, is 
not mature. The branching fractions of Y decays in- 
to BB (BB refers to both AA and X}? throughout 
the text) final states from different experiments 
and the Particle Data Group (PDG) |4] averages 
are summarized in Table [I] Obvious differences be- 
tween the different experiments are observed, and 
the uncertainties are relatively large. Hence, high- 
er precision measurements of the + decays into BB 
pairs are desirable to help in understanding the dy- 
namics of z decay. 

The angular distribution of the decays etTe~ — 
w — BB can be expressed in form 


_ By(3686)-+h 


By irish 


dN 


Teosa % 1 + acos? 0, (2) 


where 0 is the angle between the outgoing baryon 
and the beam direction in the ete~ center-of-mass 
(c.m.) system, and a is a constant, which is re- 
lated to the decay properties. The equation is de- 
rived from the general helicity formalism |], tak- 
ing into account the gluon spin, the quark distribu- 
tion amplitudes in efe — w — BB, and hadron 
helicity conservation. The e values in the decays 
J/ — BB have been calculated with pQCD to 
first-order fig. It is believed that the masses of the 
baryon and quark must be taken into consideration 
in the o calculation since a large violation of he- 
licity conservation is observed in w decays Ba) Li. 
Table [summarizes the theoretical predictions and 
experimental a values for the decays J/y — BB. 
To date, the experimental a values for the decays 
J/b — BB have poor precision (6, R! u, and 
the alpha values in the decay (3686) + BB have 
not yet been measured. It is worth noting that 
there is an indication that the o value in the de- 


cay J/y — XPS? is negative in Ref. {Li}. 
n this paper, we report precise measurements 


of the branching fractions and a values for the 
decays Y — BB, based on the data samples of 
(1310.6 + 7.0) x 10° J/w and (447.9 + 2.9) x 105 


4 


4 (3686) events collected with the BESIII detec- 
tor at the BEPCII collider. 


II. BESIII DETECTOR AND DATA SET 


The BESIII detector at the double-ring 
Beijing Electron-Positron Collider (BEPCII) is 
designed for studies of physics in the 7-charm en- 
ergy region PA. The peak luminosity of BEPCII 
is 1033 cm7? s7! at a beam current of 0.93 A. 
The BESIII detector has a geometrical acceptance 
of 93% of 4r solid angle and consists of the fol- 
lowing main components: (1) A small-celled, he- 
lium based (40% CO> and 60% C3Hg) main drift 
chamber (MDC) with 43 layers, which has an aver- 
age single-wire resolution of 135 wm, a momentum 
resolution for 1 GeV/c charged particles in a 1 T 
magnetic field of 0.5% Ï, and a specific energy 
loss (dE/dx) resolution of better than 6%. (2) An 
electromagnetic calorimeter (EMC), which consists 
of 6240 CsI (Tl) crystals arranged in a cylindrical 
shape (barrel) plus two end-caps. For 1.0 GeV pho- 
tons, the energy resolution is 2.5% (5%) in the bar- 
rel (end-caps), and the position resolution is 6 mm 
(9 mm) for the barrel (end-caps). (3) A time-of- 
flight (TOF) system, which is used for particle iden- 
tification (PID). It is composed of a barrel made of 
two layers, each consisting of 88 pieces of 5 cm thick 
and 2.4 m long plastic scintillators, as well as two 
end-caps each with 96 fan-shaped 5 cm thick plastic 
scintillators. The time resolution is 80 ps (110 ps) 
in the barrel (end-caps), providing a K/r separa- 
tion of more than 2o for momenta up to 1.0 GeV/c. 
(4) A muon chamber system, which is made of re- 
sistive plate chambers (RPCs) arranged in 9 layers 
(8 layers) in the barrel (end-caps) with ~ 2 cm po- 
sition resolution. It is incorporated into the return 
iron yoke of the superconducting magnet. 

The optimization of the event selection and the es- 
timations of the signal detection efficiency and back- 
ground are determined using Monte Carlo (MC) 
simulations. The GEANT4-based simulation 
software BOOST BA, which includes the geometric 
and material description of the BESIII detector, the 
detector response and digitization models, as well as 
the tracking of the detector running conditions and 
performance, is used to generate MC samples. The 
analysis is performed in the framework of the BESIII 
offline software system (BOSS) which takes care 
of the detector calibration, event reconstruction and 
data storage. 


+ For the J/y data sample collected in 2012, the magnetic 
field was 0.9 T. 
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TABLE I: Experimental measurements and PDG averages for the branching fractions of the decay y — BB (x10 1). 


J/b > AA (3686) > AA J/y > EE (3686) > > > 
MARKII Collab. [6] 15.8 + 0.8 + 1.9 15.8+1.6+2.5 
DM2 Collab. [7] 13.8 + 0.5 + 2.0 10.6+0.4+ 2.3 zas 
BES Collab. [8,9] 10.8 + 0.6 + 2.4 1.8 + 0.2 + 0.3 1.2 + 0.4 + 0.4 
CLEO Collab. [10] ine 3.3 + 0.3 + 0.3 sa 2.6+0.4+0.4 
BESII Collab. [11, 12] 20.3 + 0.3 + 1.5 3440.2+0.4 133404411 244+04+40.4 
BaBar Collab. [13] 19.3 +2.1 +0.5 6.41.8 +0.1 11.5 + 2.4 + 0.3 sa 
S. Dobbs et al. [14] š 9.8 = 0.1 = 0.3 s 2.3 + 0.2 + 0.2 
PDG [4] 16:1 £1.5 3:6 £0.2 12.9 + 0.9 2.3 £02 


TABLE II: Theoretical predictions and experimental 
measurements of a for J/y > BB. 


CUA AA Uy p= LEO 
0.32 0.31 [16] 
Theory 0.51 0.43 [17] 
0.72 +0.36 0.70 41.10 [6] 
Experiment 0.62 +0.22 0.22 0.31 [7] 


0.65 +0.14 —0.22 + 0.19 [L1] 


Generic inclusive MC samples, which include 
1,225 x 10° J/w and 460 x 10° (3686) events, are 
used to study the potential backgrounds. The w are 
produced via efe — ~ processes by the generator 
KKMC le, which includes the beam energy spread 
according to the measurement of BEPCII and the 
effect of initial state radiation (ISR). The known 
decay modes are generated with BesEvtGen 
according to world average branching fraction val- 
ues |M); the remaining unknown decay modes are 
simulated using the LundCharm model Bg. To de- 
termine the detection efficiencies, large y > BB 
signal MC samples are generated for each process, 
where the angular distributions of the baryons use a 
values obtained in this analysis. The A and X? par- 
ticles are simulated in the A — pr and X° > yA 
decay modes. 


II. EVENT SELECTION 


In this analysis, the four decay modes y > BB 
are studied by fully reconstructing both B and 
B, where the A(A) and ©°(S°) candidates are re- 
constructed with the pz (pr) and yA(yA) decay 
modes, respectively. Therefore, the decays y > AA 
and Y — XPE? have the final states pprtm— and 
ppr* mn yy, respectively. 

Events with at least four charged tracks with to- 
tal charge zero are selected. Each charged track 
is required to have |cos@| < 0.93, where 0 is the 


polar angle of the track. Photons are reconstruct- 
ed from isolated showers in the EMC which are at 
least 30 degrees away from the anti-proton and 10 
degrees from other charged tracks. The energy de- 
posited in the nearby TOF counters is included to 
improve the photon reconstruction efficiency and en- 
ergy resolution. Photon candidates are required to 
be within the barrel region (|cos6| < 0.8) of the 
EMC with deposited energy of at least 25 MeV, or 
within the end cap regions (0.86 < |cos6| < 0.92) 
with at least 50 MeV, where @ is the polar angle of 
the photon. In order to suppress electronic noise and 
energy deposits unrelated to the event, the timing 
information t from the EMC for the photon candi- 
date must be in coincidence with the collision event 
(0 < t < 700 ns). At least two photons are required 
in the analysis of y — ©°X° decays. 

MC studies indicate that the proton and pion 
from A decay are well separated kinematically since 
the proton carries most of the energy. A charged 
track with momentum p > 0.5 GeV/c is assumed 
to be a proton, while that with p < 0.5 GeV/c is 
assumed to be a pion. The A (A) candidate is re- 
constructed with any pz (pr) combination sat- 
isfying a secondary vertex fit a and having a de- 
cay length larger than 0.2 cm to suppress the non-A 
(non-A) decays. The decay length is the distance 
between its primary vertex and decay point to pz ` 
(prt), where the primary vertex is approximated 
by the interaction point averaged over many events. 
If more than one A (A) candidate is found, the one 
with the largest decay length is retained for further 
analysis. 


In the study of y — °° decay, a variable 


Am = 4/(May — Myo)? + (Mx, — Mgo)? is de- 


fined. All possible photon pairs are combined with 
the selected A and A candidates, and the Jı and y2 
candidates, which yield the smallest Am, are taken 
as the photons from the X? and $° decays, respec- 
tively. 

To suppress backgrounds, the AA invariant 
mass, MAA, is required to be within [3.05, 3.15], 
(2.82, 3.02], [3.63, 3.75] and [3.34, 3.61] GeV/c? for 
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the J/y > AA, J/p — X9x9, (3686) — AA 
and 7)(3686) — 5°S° decays, respectively. Here the 
mass window requirements for the individual decay 
modes are determined by MC studies. In the decays 
w — AA, the A candidate is required to have mass 
satisfying |Mp,+ — MA| < 30m, where Mj is the A 
nominal mass, and ow, is the corresponding mass 
resolution, which is 2.3 MeV/c? (4.0 MeV/c?) for 
the J/w (ap (3686)) decay. In the decays y — =°X°, 
the ©? candidate is required to have mass satisfy- 
ing |Mpr+y — Mgo| < 30m,,,, where Mso is the X° 
nominal mass, 7, 1s the corresponding mass res- 
olution, which is 4.3 MeV/c? (6.0 MeV/c?) for the 
J/w (ap(3686)). The candidates are further required 
to satisfy As050 >178° and Os050 >178.5° for the 
J/ and (3686) decays, respectively, where soso 
is the opening angle between the reconstructed X° 
and > candidates in the c.m. system. 


IV. BACKGROUND ESTIMATION 


To study the backgrounds, the same selection cri- 
teria are applied to the generic inclusive ý% MC sam- 
ples. For the decay J/h — AA, the dominant 
backgrounds are found to be J/y — AX + c.c., 
J/b — yK,Ks, and J/W — Yne with the sub- 
sequent decay ne — AA. For the decay J/y > 
LOZO, the main backgrounds are from J/y — 
AX? + c.c., J/y) — yne with the subsequent de- 
cay ne > AA, DOSO, AÈ? + cc, and J/y > 
X50520 + cc.. For (3686) — AA, the potential 
backgrounds are (3686) — ata7J/w, J/w > pp, 
4p(3686) — X20>9, and (3686) — AX + c.c.. For 
4p(3686) — X29520, the dominant backgrounds are 
from 7)(3686) > YXces,Xer > AA (J = 0,1,2) and 
(3686) => 2020,20 — An, 2° — A7’. All above 
backgrounds can be classified into two categories, 
i.e., backgrounds with or without AA in the final 
state. The former category backgrounds are expect- 
ed to produce a peak around the A/»° signal re- 
gion in the pz /pr™y invariant mass distributions 
and can be estimated, with the exclusive MC simu- 
lation samples using the decay branching fractions 
set according to the PDG 4. The additional un- 
determined decays of ne > X950_ AX? + c.c. and 
4p(3686) + AD°+<e.c. are estimated using the results 
from previous experiments for charmonium decay- 
ing to BB states (reference decays) [L], [12] [30), to 
be 1 and 0.1 times that for the decay m, — AA and 
0.1 times that for 7)(3686) — AA, respectively. The 
contributions of other decays to the peaking back- 
ground are negligible. The latter category of back- 
grounds are expected to be distributed smoothly in 
the corresponding mass distributions. 

The backgrounds from continuum QED process- 
es, i.e. ete” — BB decays, are estimated with the 


data samples taken at the c.m. energies of 3.08 GeV 
and 3.65 GeV, which have integrated luminosities of 
30 pb~! and 44 pb7! FE! Bg], respectively. By ap- 
plying the same selection criteria, no event survives 
in the selection of J/y — BB, while in the selection 
of ~)(3686) — BB, only a few events survive, and 
no obvious peak is observed in the A/D° mass re- 
gion. The contamination from the QCD continuum 
processes can be treated as non-peaking background 
when determining the signal yields. 


V. RESULTS 
A. Branching fractions 


With the above selection criteria, the distribu- 
tions of M,,-/M,,-y in a range of +8 times 
the mass resolution around the A/S? nominal 
mass in the J/7 and 7(3686) decays are shown 
in Fig. [| Clear A/=° peaks are observed with 
low background. To determine the signal yields, 
unbinned maximum likelihood fits are applied to 
Mpr-/Mpr-; with the mass of prt /prty’ restrict- 
ed to +3 times of resolution of A/=° nominal mass. 
In the fit, the A/D° signal shape is described by 
the simulated MC shape convolved with a Gaussian 
function to account for the difference in mass reso- 
lution between data and MC simulation. The peak- 
ing backgrounds are described with the shapes from 
exclusive MC simulations with fixed magnitudes ac- 
cording to the branching fractions of background 
listed in the PDG id, and the non-peaking back- 
grounds are described with second-order polynomial 
functions with free parameters in the fit. The fit re- 
sults are illustrated in Fig. [I] and the corresponding 
signal yields are summarized in Table II] 


The branching fractions are calculated using 


N. obs 


i. ST 


(3) 


where Nops is the number of signal events minus 
peaking background; < is the detection efficiency, 
which is estimated with MC simulation incorporat- 
ing the cos@ distributions obtained in this analysis 
and the scale factors to account for the difference in 
efficiency between data and MC simulation as de- 
scribed below; B; is the product of branching frac- 
tions for the intermediate states in the cascade decay 
from the PDG [4]; and Ny is the total number of 4 
events estimated by counting the inclusive hadronic 
events |18, Eg). The corresponding detection effi- 
ciencies and the resultant branching fractions are 
also summarized in Table [II] 
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TABLE III: The numbers of observed signal events Nops, the corrected detection efficiency <€, the numbers of peaking 
backgrounds Npk, the numbers of smooth backgrounds Nsm, the resultant a values for the angular distributions and 
the branching fractions B, where the errors are statistical only. 


Channel Nobs e (%) Npk Nsm a B (x10~*) 
JJ > AA. 440,675 +670 42.37 + 0.14 1,819 154+166 0.469+0.026 19.43 + 0.03 
J/w = 539559 _ 111,026 + 335 17.83 +0.06 820 131412 —0.449 + 0.020 11.64 + 0.04 
4 (3686) => AA 31,119 +187 42.83 +0.34 252 352+65 0.824+0.074 3.97+0.02 
(3686) + 530550 6,612 + 82 14.79+0.12 89 ITS 0.71 E 0.11 2.44 + 0.03 
6 3 
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FIG. 1: (color online) The M,,- distributions for the 
decays (a) J/y — AA and (b) (3686) — AA, and the 
Mpz- distributions for the decays (c) J/y — X° Ë’ and 
(d) (3686) — E° E°, where the dots with error bars are 
data, the red solid curves are the overall fit results, the 
green dashed histograms are the backgrounds estimated 
with the exclusive MC simulated samples, and the blue 
dotted line describes the remaining backgrounds. 


B. Angular distributions 


The baryon cos@ distributions in the c.m. sys- 
tem corrected by detection efficiency are shown in 
Fig. 2] and the signal yields in each of the 20 bins 
are determined with the same method as that in 
the branching fraction measurements. The detec- 
tion efficiencies in each bin are estimated with the 
signal MC samples and scaled with correction fac- 
tors to compensate for the efficiency difference be- 
tween data and MC simulation. The efficiency cor- 
rected cos@ distributions are fitted with Eq. P] with 
a least squares method, the corresponding fit results 
are shown in Fig. 2] and the resultant a values are 
summarized in Table [I] 

The correction factors used to correct for the effi- 
ciency differences between data and MC simulation 
as a function of cos are determined by studying 


FIG. 2: (color online) The distributions of efficiency 
corrected polar angle of the baryon for the decays (a) 
J/y — AA, (b) (3686) — AA, (c) J/y — XPE’, and 
(d) (3686) + °°, where the dots with error bars are 
data, and the red solid curves are the fit results. 


various control samples, where @ is the polar an- 
gle of the hyperon. The efficiency differences are 
due to differences in the efficiencies of charged par- 
ticle tracking, photon detection, and hyperon re- 
construction. For example, the efficiencies related 
with charged particle tracking and A reconstruc- 
tion are studied with a special control sample of 
w — AA events, where a A tag has been recon- 
structed. Events with two or more charged tracks, 
in which a p and z+ have been identified using parti- 
cle identification, are selected. The A tag candidate 
must satisfy a secondary vertex fit, have a decay 
length greater than 0.2 cm, and satisfy mass and 
momentum requirements. The numbers of tagged A 
events, Niag, are obtained by fitting the A peak in 
the distribution of invariant mass recoiling against 
the A tag. The numbers of A signal events, Nsig, are 
obtained by fitting the recoil mass distribution for 
events where, in addition, a A signal is reconstruct- 
ed on the recoil side, which requires two oppositely 
charged tracks that satisfy a vertex fit and have a 
decay length greater than 0.2 cm. The combined 
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efficiency of charged tracking (proton and pion) and 
A reconstruction is then Nsig/Ntag- The ratios of 
the data and MC simulation efficiencies as a func- 
tion of cos 0 are taken as the correction factors. The 
A correction factors are determined in an analogous 
way using 7 — AA events with a A tag. The over- 
all correction factor in the different cos @ bins is the 
product of the A and A correction factors. 

In an analogous way, the combined efficiency of 
photon detection and ©° reconstruction is studied 
with a control sample of y — °° events, which 
have a ©° tag and an additional A. Events are se- 
lected that have a A and A using the same criteria 
as above and at least one additional photon. The 
A and photon must have an invariant mass consis- 
tent with that of a 5°. The numbers of tagged X° 
events are obtained by fitting the X° peak in the 
distribution of mass recoiling against the X tag. 
We then search for another photon and reconstruct 
the > by requiring the invariant mass of the pho- 
ton and tagged A be consistent with the X° mass. 
The number of events with a X° signal divided by 
the number of tagged X° events is the combined ef- 
ficiency of photon detection and X reconstruction. 
The ratios of detection efficiencies in the different 
cos @ bins between data and MC simulation, deter- 
mine the correction factors. The overall correction 
factor in the different cos 0 bins is the product of the 
xO > A, and A correction factors. 


VI. SYSTEMATIC UNCERTAINTY 
A. Branching Fraction 


Systematic uncertainties in the branching frac- 
tion measurements are mainly due to the differences 
of detection efficiency and resolution between data 
and MC simulation. The sources of uncertainty re- 
lated with the detection efficiency include charged 
tracking, photon detection, and A/S? reconstruc- 
tion. The sources of uncertainty due to the resolu- 
tion difference include the MAA and Mx /Mso mass 
requirements, and the opening angle 0>oso require- 
ment in the decays y — ©°}°. Additional uncer- 
tainty sources including the model of the baryon 
polar angular distribution, the fit procedure, the 
decay branching fractions of A/X=° states and the 
total number of 7 events are also considered. All 
of systematic uncertainties are studied in detail as 
discussed in the following: 


1. As described above, the detection efficiencies 
related with the tracking, photon detection, 
and A/X° reconstruction are corrected bin- 
by-bin in cos@ to decrease the difference be- 
tween data and MC simulation. The overall 


correction factors, which are determined with 
control samples are 0.9974 + 0.0041, 0.9936 + 
0.0064, 0.980 +0.011, and 0.954 + 0.022 for the 
decays J/y — AA, J/y — X29050 (3686) > 
AA and (3686) — %°S°, respectively. To 
estimate the corresponding uncertainties, the 
correction factors are changed by +1 standard 
deviations, and the resultant changes in the 
branching fractions are taken as the system- 
atic uncertainties. 


. The uncertainties related with the MAA re- 


quirement are estimated by varying the mass 
requirement edges by +10 MeV/c?. The un- 
certainties related with the A/X° mass re- 
quirement are estimated by changing the re- 
quirement by +1 times the mass resolution. 
The uncertainties due to the requirement on 
the opening angle Os050 in the decays 7 — 
XIE? are estimated by changing the require- 
ment to be 175°. The relative changes in the 
branching fractions are individually taken as 
the systematic uncertainties. 


. MC simulations indicate that the detection ef- 


ficiencies depend on the distributions of bary- 
on polar angular cos@. In the analysis, the 
measured a values are used for the cos 0 dis- 
tributions in the MC simulation. Alternative 
MC samples are generated by changing the a 
values by +1 standard deviations and are used 
to estimate the detection efficiencies. The re- 
sultant changes in the detection efficiencies 
with respect to their nominal values are taken 
as the systematic uncertainties. 


. The sources of systematic uncertainty asso- 


ciated with the fit procedure include the fit 
range, the signal shape and the modeling of 
backgrounds. The uncertainties related with 
the fit range are estimated by changing the 
range by +1 times the mass resolution for the 
fits. The signal shapes are modeled with the 
signal MC simulated shapes convolved with a 
Gaussian function in the nominal fit. The cor- 
responding uncertainties are estimated with 
alternative fits with different signal shapes, 
i.e., a Breit-Wigner function convolved with 
a Gaussian function for A and with a Crystal 
Ball function for 5°, where the Gaussian 
function and Crystal Ball function represent 
the corresponding mass resolutions. The un- 
certainties related with the peaking back- 
grounds, which are estimated with the ex- 
clusive MC samples in the nominal fits, are 
studied by changing the branching fractions 
of the individual background, or by changing 
the branching fractions for the reference de- 
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cays which the estimated branching fractions 
for the undetermined backgrounds are based 
on, by +1 times their uncertainties from the 
PDG H. The uncertainties associated with 
the non-peaking backgrounds are estimated 
with alternative fits by replacing the second 
order polynomial function with a first. order 
polynomial function. The resultant changes 
from the above changes in the signal yields 
are taken individually as the systematic un- 
certainties. 


5. The uncertainties related with the branching 
fractions of baryon and anti-baryon decays are 
taken from the PDG H. The total numbers of 
qp events are obtained by studying the inclu- 
sive hadronic events, and their uncertainties 
are 0.6% and 0.7% for the J/w and w(3686) 
data samples {ial ud, respectively. 


The various systematic uncertainties in the 
branching fraction measurements are summarized in 
Table[IV] The total systematic uncertainties are ob- 
tained by summing the individual values in quadra- 
ture. 


TABLE IV: Systematic uncertainties in the measure- 
ment of branching fractions (%). 


IIb p (3686) 
AA. A AA 529 
Efficiency correction 0.5 0.7 1.2 2:3 
MAA requirement 0.1 0.1 0.1 0.2 
A/>° mass requirement 0.1 0.3 0.3 0.2 
Os050 requirement = 0.3 = 0.2 
Baryon polar angle 0.8 0.9 2.0 3.1 
Fit range 0.1 0.1 0.2 0.2 
Signal shape 0.1 0.3 0.1 0.2 
Peaking bkg. 0.3 0.4 0.3 1.2 
Non-peaking bkg. 0.1 0.1 0.3 0.2 
Branching fractions 1.2 1.2 1.2 1.2 
N 3y /Ny (3686) 0.6 0.6 0.7 0.7 


Total 1.7 1.9 2.8 4.3 


B. Angular Distribution 


The sources of systematic uncertainties in the 
baryon polar angular measurements include the sig- 
nal yields in different cos@ intervals and the cos 0 
fit procedure. The MC statistics and correction er- 
rors are already included in the error referred to as 
“statistical” . 


1. In the polar angular measurements, the sig- 
nal yield in a given cos@ interval is obtained 


with the same fit method as that used in the 
branching fraction measurements. The un- 
certainties of the signal yield in each cos é 
bin are mainly from the fit range, the signal 
shape and the background modeling. We in- 
dividually estimate the uncertainty of the sig- 
nal yield in each cos @ interval with the same 
methods as those used in the branching frac- 
tion measurements for the different uncertain- 
ty sources, and then repeat the cos@ fit pro- 
cedure with the changed signal yields. The 
resultant changes in the a values with respect 
to the nominal values are taken as systematic 
uncertainties. 


2. The sources of systematic uncertainty related 
to the cos @ fit procedure include the fit range 
and the number of bins in the cos@ distribu- 
tion. We repeat the fit procedures with the 
alternative fit range [—0.9, 0.9] and alternative 
number of bins (40). The resultant changes of 
qa values are taken as the systematic uncer- 
tainties. 


The individual absolute uncertainties in the po- 
lar angular distribution measurements are summa- 
rized in Table[V] The total systematic uncertainties 
are obtained by summing the individual values in 
quadrature. 


TABLE V: Absolute systematic uncertainties in the 
measurement of a. 


IIb p (3686) 

AA DEY AA ye 
Mass fit range 0.001 0.001 0.003 0.005 
Signal shape 0.001 0.002 0.001 0.003 
Peaking bkg. 0.006 0.005 0.006 0.015 
Non-peaking bkg. 0.002 0.001 0.004 0.002 
a fit range 0.001 0.003 0.007 0.019 
Number of bins 0.004 0.005 0.001 0.024 


Total 0.008 0.008 0.011 0.035 


VII. SUMMARY 


In summary, using the data samples of 1310.6 x 
10° J/ events and 447.9 x 10° (3686) events col- 
lected with the BESIII detector at the BEPCII col- 
lider, the J/y and 7(3686) decaying into AA and 
°° pairs are studied. The decay branching frac- 
tions and a values are measured, and the results 
are summarized in Table The branching frac- 
tions for J/w decays are in good agreement with the 
results of BESII and BaBar experiments, 
and those for (3686) decays are in agreement with 
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the results of CLEO [10], BESII [12] and S. Dobbs 
et al. {14] with a maximum of 2 times of standard 
deviations. The earlier experimental results 

have significant differences with those of this analy- 
sis. The precisions of our branching fraction results 
are much improved than those of previous experi- 
ments listed in Table [| The a values in the decays 
w(3686) > AA and (3686) — ©°S° are measured 
for the first time, while those of J/y — AA and 
J/y — XPE? decays are of much improved precision 
compared to previous measurements. It is worth 
noting that the a value in the decay J/y — X09050 
is negative, which confirms the results in Ref. H. 


TABLE VI: Results for measured o values and branch- 
ing fractions ËB in this analysis. The first uncertainties 
are statistical, and the second are systematic. 


Channel a B (x107*) 
Jp > AK 0.469 + 0.026 £ 0.008 19.43 £0.03 £0.33 
J/h — ELS? 0.449 + 0.020 + 0.008 11.64 + 0.04 + 0.23 


w(3686) > AA 0.82 + 0.08 + 0.02 3.97 + 0.02 + 0.12 


(3686) — XPE? 0.71 40.1140.04 2.44 + 0.03 + 0.11 


To test the “12% rule”, we also obtain the Q val- 


ues to be SED S = (20.43 + 0.11 + 0.58)% 
and 800) — (20.96--0.27--0.92)%, where 


the common systematic uncertainties between .J/*) 
and ~(3686) decays are cancelled. The Q values are 
of high precision, and differ from the expectation 
from pQCD by more than 3 standard deviations. 


VIII. ACKNOWLEDGMENTS 


The BESIII collaboration thanks the staff of 
BEPCII and the IHEP computing center for their 


1] P. Kessler, Nucl. Phys. B 15, 253-266 (1970); 
S. J. Brodsky and G. P. Lepage, Phys. Rev. D 24, 
2848 (1981). 

2| T. Appelquist and H. D. Politzer, Phys. Rev. Lett. 
34, 43 (1975); A. De Rujula and S. L. Glashow, 
Phys. Rev. Lett. 34, 46 (1975); W. S. Hou, Phys. 
Rev. D 55, 6952 (1997). 

3] M. E. B. Franklin et al. (MARKII Collaboration) 
Phys. Rev. Lett. 51, 963 (1983). 

4] K. A. Olive et al. (Particle Data Group), Chin. 
Phys. C 38, 090001 (2014). 

5| Y. F. Gu and X. H. Li, Phys. Rev. D 63, 114019 
(2001); X. H. Mo, C. Z. Yuan and P. Wang, High 
Energy Physics and Nuclear Physics 31, 686 (2007); 
N. Brambilla et al. (Quarkonium Working Group), 


10 


strong support. This work is supported in part by 
National Key Basic Research Program of China un- 
der Contract Nos. 2009CB825200, 2015CB856700; 
National Natural Science Foundation of China 
(NSFC) under Contracts Nos. 10905034, 10935007, 
11125525, 11235011, 11322544, 11335008, 11425524; 
the Chinese Academy of Sciences (CAS) Large- 
Scale Scientific Facility Program; the CAS Center 
for Excellence in Particle Physics (CCEPP); the 
Collaborative Innovation Center for Particles and 
Interactions (CICPI); Joint Large-Scale Scientific 
Facility Funds of the NSFC and CAS under 
Contracts Nos. 11179007, U1232106, U1232201, 
U1332201; Natural Science Foundation of Shandong 
Province under Contract No. ZR2009AQ002; CAS 
under Contracts Nos. KJCX2-YW-N29, KJCX2- 
YW-N45; 100 Talents Program of CAS; National 
1000 Talents Program of China; INPAC and 
Shanghai Key Laboratory for Particle Physics and 
Cosmology; German Research Foundation DFG un- 
der Contract No. Collaborative Research Center 
CRC-1044, FOR-2359; Istituto Nazionale di Fisica 
Nucleare, Italy; Joint Funds of the National 
Science Foundation of China under Contract No. 
U1232107; Ministry of Development of Turkey un- 
der Contract No. DPT2006K-120470; Russian 
Foundation for Basic Research under Contract 
No. 14-07-91152; The Swedish Resarch Council; 
U. S. Department of Energy under Contracts 
Nos. DE-FG02-04ER41291, DE-FG02-05ER41374, 
de-sc0012069, DESC0010118; U.S. National Science 
Foundation; University of Groningen (RuG) and 
the Helmholtzzentrum fuer Schwerionenforschung 
GmbH (GSI), Darmstadt; WCU Program of 
National Research Foundation of Korea under 
Contract No. R32-2008-000-10155-0. 


Eur. Phys. J. C 71, 1534 (2011); Q. Wang, G. Li 
and Q. Zhao, Phys. Rev. D 85, 074015 (2012). 
[6] M. W. Eaton et al. (MARKII Collaboration), Phys. 
Rev. D 29, 804 (1984). 
[7] D. Pallin et al. (DM2 Collaboration), Nucl. Phys. 
B 292, 653 (1987). 
[8] J. Z. Bai et al. (BES Collaboration), Phys. Lett. B 
424, 213 (1998). 
[9] J. Z. Bai et al. (BES Collaboration), Phys. Rev. D 
63, 032002 (2001). 
[10] T. K. Pedlar et al. (CLEO Collaboration), Phys. 
Rev. D 72, 051108 (2005). 
[11] M. Ablikim et al. (BESI Collaboration), Phys. 
Lett. B 632, 181 (2006). 
[12] M. Ablikim et al. (BESII Collaboration), Phys. 


201711.02285v1 


chinaXiv 


13 


14 


15 


16 


17 
18 


[19] 


[20] 


Lett. B 648, 149 (2007). 

B. Aubert et al. (BaBar Collaboration), Phys. Rev. 
D 76, 092006 (2007). 

S. Dobbs et al., Phys. Lett. B 739, 90 (2014). 

M. Ablikim et al. (BESIII Collaboration), Phys. 
Rev. D 93, 072003 (2016); M. Ablikim et al. (BESIII 
Collaboration), arXiv:1612.08664. 

M. Claudson, S. L. Glashow and M. B. Wise, Phys. 
Rev. D 25, 1345 (1982). 

C. Carimalo, Int. J. Mod. Phys. A 2, 249 (1987). 
M. Ablikim et al. (BESIII Collaboration), Chin. 
Phys. C 36, 915 (2012); M. Ablikim et al. (BESIII 
Collaboration), arXiv:1607.00738, (to be published 
on Chin. Phys. C). 

M. Ablikim et al. (BESIII Collaboration), Chin. 
Phys. C 37, 063001 (2013); The total number of 
~)(3686) events taken at 2009 and 2012 is obtained 
based on the same method. The preliminary num- 
ber is determined to be 447.9 x 10° with uncertain- 
ties of 0.7%. 

M. Ablikim et al. (BESIII Collaboration), Nucl. 
Instrum. Meth. A 614, 345 (2010). 


21 
22 


23 


24 
25 


26 


27 


28 


29 
30 


31 


11 


C. Zhang, Sci. China G 53, 2084 (2010). 

D. M. Asner et al., Int. J. Mod. Phys. A 24, 499 
(2009). 

S. Agostinelli et al. (GEANT4 Collaboration), Nucl. 
Instrum. Meth. A 506, 250 (2003). 

Z. Y. Deng et al., HEP & NP 30, 371 (2006). 

W. D. Li, Y. J. Mao and Y. F. Wang, Int. J. Mod. 
Phys. A 2481, 9 (2009). 

S. Jadach, B. F. L. Ward and Z. Was, Comp. Phys. 
Commu. 130, 260 (2000); Phys. Rev. D 63, 113009 
(2001). 

D. J. Lange, Nucl. Instrum. Meth. A 462, 152 
(2001); R. G. Ping, Chin. Phys. C 32, 599 (2008). 
J. C. Chen, G. S. Huang, X. R. Qi, D. H. Zhang 
and Y. S. Zhu, Phys. Rev. D 62, 034003 (2000). 
M. Xu et al., Chin. Phys. C 33, 428 (2009). 

M. Ablikim et al. (BESIII Collaboration), Phys. 
Rev. D86, 032014 (2012). 

J. Cheng, Z. Wang, L. Lebanowski, G. Lin and 
S. Chen, Nucl. Instrum. Meth. A 827, 165 (2016). 


